Evidence for a conformal phase in SU(N) gauge theories 
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We discuss the existence of a conformal phase in SU(N) gauge theories in four dimensions. Simu- 
lations are carried out with three colors and twelve flavors of dynamical staggered fermions in the 
fundamental representation. The analysis of the chiral order parameter and the mass spectrum of 
the theory indicates the restoration of chiral symmetry at zero temperature and the presence of a 
Coulomb phase, depicting a scenario compatible with the existence of an infrared stable fixed point 
at nonzero coupling. Our analysis thus supports the conclusion that the onset of the conformal 
window for QCD-like theories is smaller than Nf = 12, before the loss of asymptotic freedom at 
sixteen and a half flavors. 



With the imminent activity of the LHC experiments 
and the quest for a theory describing fundamental forces 
beyond the electroweak symmetry breaking scale, re- 
newed interest has arisen in the most elusive aspects 
of gauge theories. In particular, the possibility of an 
emergent quasi-conformal symmetry in theories with 
fermionic content has attained a strong experimental ap- 
peal. There are multiple reasons for pursuing this search. 
Resolving conformal behavior would complete our under- 
standing of the phase diagram of gauge theories by vary- 
ing temperature and number of flavors, as sketched in 
Fig. [T] It sheds light on how the low temperature and 
large flavor number quasi-conformal phase may be con- 
nected to the high temperature and low flavor number 
quark-gluon plasma phase. It is essential for theoretically 
establishing or excluding walking technicolor-type the- 
ories and more generally strongly interacting dynamics 
above the electroweak symmetry breaking scale. Finally, 
elucidating the way conformal symmetry or its remnants 
drive the dynamics of particle interactions with or with- 
out supersymmetry contributes to clarifying the possi- 
ble connection of field theory to string theory that the 
AdS/CFT correspondence seems to imply. 

In the early eighties our understanding of the pertur- 
bative behavior of non Abelian gauge theories was en- 
riched by two seminal papers [TJ [2] . It was noticed that 
a second zero of the two-loop beta function of an SU(3) 
gauge theory with Nf massless fermions in the funda- 
mental representation appears for Nf > 8.05, at g* 2 ^ 0, 
before the loss of asymptotic freedom at Nf = 16 1. This 
fact implies, at least perturbatively, the appearance of an 
infrared fixed point (IRFP). The fixed point moves closer 
to zero coupling as the number of flavors approaches Nf. 
The dynamics of chiral symmetry have led to the discov- 
ery of the conformal window in QCD-like theories [3111]: 
chiral symmetry breaking can only occur below a critical 
number of flavors Nf . Between Nf and Nf the confor- 
mal window opens up. Finding the actual value of the 
critical number of flavors Nf at which chiral symmetry 
breaking takes place is a non-perturbative problem for 
which the lattice approach is the method of choice. Re- 
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FIG. 1: A projected view of the phase diagram of QCD-like 
theories in the temperature (T), flavor number (Nf) and bare 
coupling (g) space. In the T-Nf plane, the critical line is a 
phase boundary between the chirally broken hadronic phase 
and the chirally symmetric quark gluon plasma, the zero tem- 
perature end point of which is the onset of the conformal win- 
dow. The zero temperature projected plane is inspired by the 
scenario in Refs. [3j [4] , see Fig. [2] 

cent studies have provided evidence that Nf = 8 lies 
within the hadronic phase of QCD EJ [9] . A re- 
cent study of the SU(3) running coupling [SI IH] by use 
of the lattice Schrodinger functional has concluded that 
Nf = 12 should already be in the conformal window. 
Other numerical studies, however, challenged this con- 
clusion [lOj HH [12] . This is hardly surprising, given that 
Nf = 12 should be very close to the critical number of 
flavors [T31 HH E3 HH] , making a numerical study partic- 
ularly delicate. 

The current strategy, complementary to that of 
Refs. [5J [5], is inspired by the physics of phase transi- 
tions; it allows for the exploration of multiple aspects of 
the theory in different regimes and regions of the phase 
diagram, in order to probe the existence and properties 
of an IRFP inside and outside its basin of attraction. 

A lattice strategy The strategy of this study has 
received heuristic guidance by the scenario depicted in 
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FIG. 2: (color online) Phase diagram of an 5)7(3) gauge the- 
ory with fundamental fermions in the Nf-g plane after Ref. [4] . 
Theories for Nf < Nf are QCD-like in the continuum, while 
for Nf < Nf < Nf develop a conformal phase. S and A re- 
fer to chirally symmetric and asymmetric, respectively. The 
dashed(green) line qualitatively indicates the location of the 
Banks-Zaks IRFP [2]. The dot-dashed (red) line indicates a 
lattice bulk transition, which has been observed at Nf = 12 
and Nf = 16. It stops at some small Nf value. The line 
at Nf — Nf represents the conformal phase transition [3] 2] , 
which is absent in the original Banks-Zaks scenario. The beta 
function on the conformal side is also sketched. 



Refs. [3j 2] and sketched in Fig. [2j We will implicitly as- 
sume the validity of that scenario. Following Fig. [2] at a 
given Nf > and increasing the coupling from g = 0, 
one crosses the conformal line, location of the IRFP's, go- 
ing from a chirally symmetric (S) and asymptotically free 
phase (quasi-conformal phase) to a symmetric, but not 
asymptotically free one (Coulomb or QED-like phase) . A 
phase transition need not be associated with the line of 
IRFP's, differently from what was originally speculated 
in Ref. [2 . At even larger couplings, a transition to a 
strongly coupled chirally asymmetric (A) phase will al- 
ways occur in the lattice regularized theory. The latter 
is referred to as a bulk phase transition. In the symmet- 
ric phases at nonzero coupling the conformal symmetry is 
still broken by ordinary perturbative contributions. They 
generate the running of the coupling constant which is 
different on the two sides of the symmetric phase. See 
Ref. jl] for a detailed discussion of this point. 

A theory in the hadronic phase, Nf < NJ-, has a ther- 
mal phase transition in the continuum from a low tem- 
perature chirally broken phase to a high temperature chi- 
rally symmetric - quark gluon plasma - phase. Thus, as 
argued in Ref. [6], the observation of a thermal tran- 
sition in the continuum limit is incompatible with the 
existence of a conformal fixed point, see Fig. [T] It is also 
clear from Fig. [2] that the presence of a Coulomb phase 
next to the bulk transition at weaker coupling is also a 
distinguishing feature of the conformal window. Here, 
the non-perturbative beta function should be positive, 
implying a weakening of the effective coupling over in- 
creasing distances. The appearance of such a region is, 
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FIG. 3: (color online) The bulk transition in the chiral con- 
densate for am = 0.05 on lattices of 16 3 x 8 (circles), and 16 4 
(crosses) as a function of the bare lattice coupling gL- Data 
are shown in the range 6/g'l = 2.5 to 4.7. The location of 
the transition is identical, while the curves describe physics 
at temperatures differing by a factor of two. Simulation errors 
are within symbol size. 



in principle, a sufficient condition for the existence of an 
IRFP, since the perturbative beta function of 5(7(3) with 
Nf < 16 1 in the extreme weak coupling regime is known 
to be negative. Note, however, that the beta function is 
not universal away from fixed points with diverging corre- 
lation lengths and one can therefore not exclude a priori 
the appearance of spurious fixed points at intermediate 
values of the coupling constant [17] . 

The evidence presented here thus consists of a few com- 
ponents. First, it will be demonstrated that the location 
of the transition from the chirally symmetric to the bro- 
ken phase is not sensitive to the physical temperature and 
is therefore compatible with a bulk nature. Subsequently, 
we will present a detailed study of the mass dependence 
of the chiral condensate on the weak coupling side of the 
bulk transition, which clearly favors exact chiral symme- 
try. Finally, the behavior of the mass spectrum close 
to the bulk transition will be studied, and found to be 
compatible with a positive beta function, similarly to the 
observations of Ref. [17] for Nt = 16. These results are 
consistent with the scenario for conformality of Fig. [2] 

Results We performed simulations of an 5(7(3) gauge 
theory with twelve flavors of staggered fermions in the 
fundamental representation. The simulations were done 
using a tree level Symanzik improved gauge action to sup- 
press lattice artifacts, and Kogut-Susskind (staggered) 
fermions with the Naik improvement scheme, that effec- 
tively extends the Symanzik improvement to the matter 
content. The order parameter of the chiral phase transi- 
tion for a QCD-like theory is the chiral condensate (iJip)- 
We underscore that staggered fermions have a remnant 
of exact chiral symmetry which allows a precise definition 
of the order parameter also on a coarse lattice. Fig. ^dis- 
plays a sudden variation of the chiral order parameter as 
a function of the bare lattice coupling constant gL, and 
for two values of the temporal lattice extent Nt differing 
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FIG. 4: (color online) Deviations of the predicted value for 
an extrapolation to the chiral limit of the chiral condensate 
for 6/gx, = 3.9 (circles) and &/g\ = 4.0 (triangles), fitting to 
a linear form Am + B (red) - with m the bare mass in lattice 
units - and a power law Am 6 (blue), rescaled by the standard 
deviation of each measurement. The data and corresponding 
fits are displayed in the inset on linear scales. 



by a factor of two. The temperature of the system is 
related to the lattice temporal extent as T — l/a(gL)N t , 
with a(<7z,) the lattice spacing for a given lattice cou- 
pling. From Fig. [3] one infers that the transition - or 
rapid crossover - happens at identical values of the crit- 
ical coupling g c L — 1.35(3), thus implying they occur at 
vastly different physical temperatures. One concludes 
that the observed phase transition is driven purely by the 
bare coupling constant itself and is therefore of bulk na- 
ture. Further information on this phase transition might 
shed light on the occurrence of a conjectured ultraviolet 
fixed point at strong coupling in the continuum theory 

To discriminate between a weakly broken and an ex- 
act chiral symmetry on the weak coupling side of the 
bulk transition, we performed a high statistics study of 
the mass dependence of the chiral condensate for two 
selected couplings. Measurements were performed on 
three different volumes for each mass, up to 32 4 , and 
the difference between the largest two volumes found to 
be smaller than both the difference between the smallest 
volumes and the statistical uncertainty in all measure- 
ments. Data displayed in Fig. [4] all at volume 24 4 , can 
be considered as infinite volume measurements within 
their errors. We fitted the chiral condensate data to 
< tj)ip >— Am+ < ijjip >o, appropriate for a chirally 
broken (< iptp >ot^ ) or symmetric phase, away from 
criticality and in the absence of an IRFP, and to a power 
law < {/ji/j >— Am s . The latter should be expected in a 
symmetric phase and in the vicinity of a fixed point with 
anomalous dimensions. For instance, this behavior was 
observed for lattice QED in Ref. [19] . The pattern of the 
deviations in the linear fits indicates a significant curva- 
ture, which is reflected in the quality of the fit. The x 2 
per d.o.f. was 0.57 (power) as compared to 3.21 (linear) 
for <7l = 1.240 (6/cr^ = 4.0), and 0.86 as compared to 



3.81 for g L = 1.224 (6/g 2 L = 3.9), with best fitted expo- 
nent 8 = 0.965(1) for both. A power law fit with a free 
intercept always produced a best fitted value consistent 
with zero. To complete our study of the symmetric phase, 
following Ref. [T7], we measured the mass spectrum and 
used the results to determine the lines of constant physics 
in the two dimensional parameter space gL and am, the 
bare quark mass of degenerate fermions. Along these 
lines the coupling and masses are all functions of the lat- 
tice spacing a. Since all dimensionful quantities measured 
on the lattice will be expressed in terms of the lattice 
spacing and will therefore vary with even if they do 
not physically, a dimensionless quantity has to be taken 
as a reference. A convenient choice is the ratio of the 
ground state excitations in the pseudoscalar and vector 
channels, with slight abuse of nomenclature from QCD 



referred to as the tt and p masses. In Fig. 5(a) we report 



on the measured values of and m p as a function of 
the bare coupling and the bare mass. Simulations were 
done on 16 3 x 24 volumes, while a set of measurements 
at larger volumes showed that finite volume effects were 
under control. 

In Fig. |5(b)| we show the ratios of measured pseu- 
doscalar and vector masses, for a fixed coupling and as 
a function of the bare quark mass. We have superim- 
posed the ratios of the best fits to the raw mass data, 
confirming the good quality of the interpolations derived 
in Fig. 5(a) and later used to produce Fig. [6] It is imme- 
diately evident that, in order to keep the ratio constant, 
we should simultaneously decrease the lattice coupling 
gL and increase the bare mass. These results already 
indicate that the lattice spacing increases while decreas- 
ing the coupling. This is the same behavior as observed 
for Nf — 16, and of the pion to sigma ratio in QED. It 
is also expected of a one parameter beta function with 
a positive sign. To refine the analysis, and express the 
result in terms of a physical observable, we proceed as 
follows. Given the ratio m^jm p at reference values of 
gL and am, one can determine a value aim' at coupling 
g' L in the surroundings of gL that reproduces the same 
ratio and thus lies on the same line of constant physics. 
This is implemented by fitting the measured values of 
both masses to a parametrization, as shown in Fig. 5(a)| 
then determining the isolines from the ratio of the two 
parametrizations. The physical pseudoscalar mass be- 
ing constant along each of these lines, the ratio of mea- 
sured values am^/a'm^ directly determines the ratio of 
the lattice spacings a/a 1 . If a decrease in gL is associ- 
ated with an increase in the lattice spacing, the sign of 
the beta function is positive, i.e. that of the Coulomb 
phase in Fig. [2] Generalizing, along the lines of constant 
physics the slope of the line of measured values of the 
pseudoscalar mass is a direct measure of the sign of the 
beta function. Fig. [6] provides evidence for the Coulomb 
phase, with a positive sign of the beta function, in full 
agreement with the more naive discussion of Fig. |5(b)| 
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(a)Interpolations of meson masses. 
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FIG. 5: (color online) (a) Measurements of the pseudoscalar (blue) 
values of the bare quark mass, from bottom to top am — 0.04, 



and vector (red) masses versus lattice coupling at several 
0.05, 0.06 and 0.07. Lines displayed represent a global 



parametrization, with a mixed 0(rn) polynomial quark mass dependence and 0(/3 ) polynomial dependence, with lattice 
parameter (3 — Q/g\, and producing a reduced \ 2 P er d.o.f. just over unity for both channels. Errors include fitting systematics 
from combining several methods, (b) The measured n to p mass ratio as a function of the b are m ass and decreasing coupling 
g L , bottom to top G/g\ = 3.5 to 4. The superimposed lines are ratios of the best fits in Fig. 5(a) 
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FIG. 6: Pseudoscalar mass along lines of constant physics. 
The physical pion mass is identical along lines of constant 
physics, such that a decreasing pion mass in lattice units for 
increasing gi, implies an increase of the lattice spacing for 
weaker couplings, signature of the Coulomb phase. The pseu- 
doscalar mass along lines of constant physics was constructed 
by interpolating along the isolines of the ratio of the separate 
interpolations of pseudoscalar and vector masses. Different 
polynomial interpolations produced compatible results, and 
in agreement with a non-interpolated analysis of raw data. 
Labels give the value of the ratio m^/mp along the isolines. 

To summarize, for an SU(3) gauge theory with twelve 
flavors the lattice phase transition from a chirally broken 
to a symmetric phase is of a non-thermal nature, as is 
expected for a theory within the conformal window. A 
study of the spectrum in the symmetric phase close to 
the transition shows a positive sign of the beta function, 
incompatible with asymptotic freedom and implying the 
emergence of a Coulomb phase. Assuming the validity of 
the scenario in Refs. [31 13] and Fig. [2] the observation of 
a Coulomb phase points to the existence of an infrared 
stable fixed point for the theory with twelve flavors and 
of a conformal window. 
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